HIGH-FREQUENCY GRAVITATIONAL WAVES FROM 
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We investigated the stability of the spherically symmetric non-abelian (Bartnik-McKinnon) black 
hole solution of the SU(2) Einstein- Yang-Mills system using the multiple-scale analysis. It is found, 
in contrast with the vacuum situation, that a spherically symmetric oscillatory perturbation to 
second order cannot be constructed. The singular behavior of the gravitational waves is probably 
induced by the coupling of the gravitational waves to the Yang-Mills waves. 
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■- I. INTRODUCTION 

ON 

On 

^\ ' The discovery of the non-abelian particle-like and non-Reissner-Nordstrom (hairy) black hole 

solutions of the SU(2) Einstein- Yang-Mills (EYM) theory, first studied by Bartnik and McKinnon 

(BMK) jjj], initiated a throughout instability investigation by several authors |^|-|^|. It turns out that 

under spherically symmetric perturbations the BM solution is unstable. The nth BMK solution has 

in fact 2n unstable modes, comparable with the flat spacetime "sphaleron" solution. In the static 

case, it was recently conjectured Q that there are hairy black holes in the SU(N) EYM-theories with 

l/") . topological instabilities. In order to analyze the stability of the solutions, one usually linearizes the 

field equations, or one expands the field variables in a physically unclear small parameter H^j. We 

conjecture that the conventional linear analysis is in fact inadequate to be applied to the situation 

where a singularity is formed. One better can apply the so-called multiple-scale or two-timing 

method, developed by Choquet-Bruhat This method is particularly useful for constructing 

uniformly valid approximations to solutions of perturbation problems JhJ. One of the advantages 

, is that one can keep track of the different orders of approximations. The method originates from 

the Wentzel-Kramers-Brillouin(WKB)-method in order to construct approximate solutions of the 

HN ' non-linear Schrodinger equation. It can also be applied when there arise secular terms, as in the 

JZL. ' .2 

■ Mathieu equation + (a + 2e cos t)y = 0. As a nice application of the multiple-scale-method, one 

""^ | finds the boundaries between the regions in the (a, e) plane for which all solutions are stable and the 

^J^i regions in which there are unstable solutions. Quite recently it was found that the method can be 

successful for gravity with Gauss-Bonnet terms. There is another application, namely the threshold 

of black hole formation of the choptuon [pd[-|l4|. Numerically it is found that there is a critical 

parameter whose value separates solutions containing black holes from those which do not. For 

a critical value one observes self-similarity: it produces itself (echoes) on progressively finer scales 

(Choptuik scaling). It is evident that the collapsing ball of field energy will produce gravitational 

waves, which will be coupled to the YM-field perturbations. Due to the accumulation of echoes, the 

curvature diverges. The multiple-scale method is suitable for handling this critical situation. 

Let us consider a manifold M with two different scales, i.e., a mapping from A4 x TZ into the 

space of metrics on M: 

(x,f)=>9(x,0, xe.M,£eft. (i) 

We set £ = cjII(x), with IT a phase function on M of dimension of length and u a large parameter 
of dimension {length)' 1 . 

The multiple-scale method assumes that the metric g^ v and the YM-potentials A a ^ can be written 

as 

gy.v = g^u + — h^x";^) + \k^{x a ;£) + (2) 
A«=A; + -B^(x'i + \c«{x°; + (3) 



where £ = ujH{x") and II a phase function. The parameter uj measures the ratio of the fast scale to 

an 



the slow one. The rapid variation only occur in the direction of the vector l a = For a function 



^(x"\C) one has 



1 



where 8^ = JrlHf fixed an( i * = ^\x a fixed- We consider here the case where the magnitude of 
the perturbation of the metric with respect to the background is of order (for different 
possibilities, such as for example with the leading term in the metric expansion of order ui~ 2 , see 
Taub |^j). Further, we will consider here the hypersurfaces Yl(x a ) — est as wave fronts for the 
background metric g, so (Eikonal equation) 

lalpg™ 13 = 0. 

Substituting the expansions of the field variables into the equations and collecting terms of equal 
orders of ui, one obtains propagation equations for C^, and and 'back-reaction' equations 
for g^ v and A^. It will be clear from the propagation equation that there will be a coupling between 
the high-frequency gravitational field and the high-frequency behavior of when the singularity 
will be approached. 



II. THE EYM FIELD EQUATIONS IN THE MULTIPLE-SCALE FORMULATION 



Consider the Lagrangian of the SU(2) EYM system 

TZ 1 



5' : 



d xy^g 



spa spf_Li>a 



with the YM field strength 



f a — ft A a — rl A a 4- nf c A" A c 



g the gauge coupling constant, G Newton's constant, A^ the gauge potential, and 71 the curvature 
scalar. The field equations then become 



Guv — — &7vGTuv , 



V^ va = 0, 

with T the energy-momentum tensor 



and D the gauge- covariant derivative, T> a J-^„ = V a ^-"^ + ge ahc A^T^u- Substituting the expansions 
into the YM-equation Eq. 9, we obtain for the order ui equation 



9 



I filaBjy It/lot B „ 



= 0. 



For the order adequation we obtain 



V7 -7~Q I *y>A A s-p< 



+ge abc (l a (Alm + AlB*) + A" a (j^„ + l»K - - b? a l a {l»Bl ~ l v B$ = 0, 

with = \g' yX {l l _ i h l/cr - Ivhua - Ivh^u)- 

Substituting the expansions into the Einstein equations Eq. 8, we obtain for the order u equation 

p(-l) _ ; y° _ l ye _ n 
li '/Liy — I/it ± fi G la ± /_iiy — U. 

For the u;°-equation we obtain 



■ 1 — —era — TjS 

+j9^g g ^ 



with i? M „ the background Ricci tensor and R$ an expression in k, hh, h 2 and h (see Choquet-Bruhat 
|7[). We can simplify the equations considerable. If we use Eq. 5, then Eq. 11 becomes 



so for periodic we have 



FBI = 0, 



FBI = 0. 



Using Eq. 13 and Eq. 15 we obtain from Eq. 12 



abc—f-Lot 



0, 



with V the covariant derivative with respect to the background metric <? M „. Integrating this equation 
with respect to £ yields 



£> M J% = 0. 



Substituting back this equation into Eq. 12 we obtain the propagation equation for the YM-field 

V"(i„Bj - l„BZ) + l»(V^BZ - V V B«) - g» a r«A + ^X-Ba 
-IXC" + ge ahc t a (2l a AlB c v - LA b a B^) = 0. 



Multiplying this propagation equation Eq. 19 with B va , we obtain the 'conservation'-equation 

V a {rKB ua ) = Fh x B va ^ x - 2ge abc l" B c v B ua A*, 
so B„B va is not conserved, unless the right-hand side is zero. Multiplying Eq. 19 with V , we obtain 

i"{i x ht-rhl)^ x = o. 

For the choice l v h vii = and hence from Eq. 13 h — 0, the equations simplify again. 
Integrating Eq. 14 with respect to £, we obtain 

5 M „ = —Ink /" {KK - - —^—Iftlv J B aa B a a dt - 8ttG%„, 

with r the period of high-frequency components. In the right-hand side we have the back-reaction 
term due to the high-frequency gravitational perturbation, the high-frequency YM-perturbations 
and the background energy- momentum term T. Substitution back Eq. 22 into Eq. 14, we obtain 
the perturbation equations for h^v 



i? (0) - I I 



1 

2 ! 



8nG 



(l IJ .l l/ B Xa B x 



+^ x (l v B Xa - l x BZ)+^ x (l^B Xa - l x B a J + y^F aP (l a B^ a - fB aa " 



III. APPLICATION TO THE SPHERICAL SYMMETRIC SPACETIME 



Let us consider the spherical symmetric spacetime 

1 



ds = -m{t,r)N(t,r) dt + 



dr + r (dO + sin Odtp 



We have for the tetrad component: i M = (- 
SU(2) YM field in the 'polar' gauge as § 



m(t, r) 

N,±,Q, 0) and Z M = 1, 0, 0). We parameterize the 



I = A^dx^ = A*T a dx>* = AiT 3 dt + +A 2 T 3 dr + (Wm + W 2 T 2 )d6> + (cot <9r 3 + Wir 2 - WaTi) sin 0d<^, 
with r a the spherical projections of the SU(2) generators. In components 

A t — 0, 0, 0, + — Oj + + 



A e = 

A<o — 



0,0,b + -bi + -U2 + .. 

UJ UJ 2 

-{w H — wi H 5-w 2 + ...), -(v(u) H — ui H + • ••), 

UJ UJ 2 UJ UJ 2 

'- 1 1 , 1 1 , . „ 

v H — ui H «2 + •■, (— (w H — toi H 7«i2 + ...) sin6», — cos( 

UJ UJ UJ UJ 2 



where the perturbations depend on t, r, 6, ip and £. From Eq. 16 we obtain the condition bi — — t^jt- 
From Eq. 13 we obtain some restrictions on h uv : hni = — „ h ° » r — i\mNh\\ , h\2 = and 

n ui 2mN 2 ' mJV 

fti3 = — ;^v- We can now evaluate the propagation equations and back-reaction equations. From 
the YM propagation equation Eq. 19 we obtain the set (g = — 1) 



w 2 + v 2 = 0, 



6 = 



miV' 

(9tWi + mNdrWi — 0, 



ftm — m d r N — 0, 



^ • • • ,dtm dtN . , t2 ;- 

m iy 

. , Q _ . Ara - ad t m ad t N (h Q0 -m 2 N 2 h±±) 
+(d t a + mNd r a — — )(- - — — = 



N 



2mN 2 



The third equation in Eq. 27 is consistent with the result that the divergence of the null- 
congruence of F must be equals 2, an d the fact that d r tH = dtrTl. The back-reaction equation 
Eq. 18 for the background YM field yields (from now on we set v — w = 0) 

«2i c.j.dtN _d t N _ 

* ' ~/v" ra_ /v — = 

For the propagation equations Eq. 23 of the high-frequency gravitational perturbations h and k, we 
obtain 



dtfi22 + mNd r h 2 2 



mNh 22 



0, 



d t h,23 + mNd r h,2z — mNfl23 — q ; 



;■ 1 ;• dth 03 , 9r/l03 

mN mN 2 N 



d t m d t N 2 \- 

~ H — H ~ I ^03 



+ 



m 2 N 2 mN 3 rN 
d v h 22 deh,23 /123 cot 6 IGtxGcosO 



d v hoo mdiphii 
2mN 2 + 2 



% 2 N 3 



(mNdta — aNdtm — amd t N + m N d r a), 



i j. -Li ( 2 i dt7Tl • dtN 

12 + Mn 02 _ \rN m 2 N 2 + mN 3 
1 



• 167rGcot0 1 ■ 2cot6»- 
no2 7, rrae/122 ; — h.22 

., . 2-8^23 + ^dehu - ^-d r h 2 Ld t h 2 - - 1 . r2 dghoo, 

r 2 sm 2 2 N mN 2 2mN 2 



k 2 2 + 



k33 



2 2 
2m rftn — 327rG(u 1 + H — ^h22h,22 + 



sin a_ 

h 2 2 ^23 
"I T ^ ,-2/1 ■ 2 n U <P n 03 

r 2 r 2 sm 2 N sm 2 6 
and 



■■ 2 • 2cotg - 
2 „ o "23^23 — —deno2 ^ — n-02 



^03 + 

r 2 r 



J2 . ^23 
/l 22 + ~. o" 



sm-^r 2 AT AT 

327fG /" , . 2 . .2s 

H ^— / (»i +vi)<^, 



2 •■ 1 

H AT kQ1 H 2TF2 

mN m 2 N 2 



koo 



( 



+ (^779 t m + d r m) hn + -^d t hu + md r h 
y mN J N 

From the back-reaction equations Eq. 22 we obtain 



2 3 1 \ • 

d t N + ,,„ 9tm+ * TO d r m\h 00 



m 2 iV 4 



m 3 iV 3 
11 „ ,„ dthoo 



2 N 2 



m 2 N 3 



mN 2 



d r hoo- 



d t m = 



N 
2rr 3 



7/4, sm 2 e + h 2 2 3 



167rGAf 



sin 2 6 + vl sin 2 6» + 5 cos 2 1 



which represents the effect on the decrease of the metric ("mass") component M due to the high- 
frequency perturbations, where M is the well known mass parameter in the Schwarzschild parame- 
terization m = (1 — — ). 



IV. CONCLUSIONS 



Let us try to find a bounded solution to second order. Following the vacuum situation of Choquet- 
Bruhat jjj], we consider the simplified situation, where we take /102 = ^03 = ko2 = ko3 = and h^v 
independent of (p. Further, 

&33 = — sin 2 6k 2 2 

h 2 2 = rH.22 (r, t) cos £ sin 8 

h,23 = rH,2z (r, t) sin £ sin 8 cos 8. (38) 

From Eq. 31, Eq. 32 and Eq. 35 we could find a bounded (not spherically symmetric) solution 
for the other components of /i M „. However, from Eq. 34 we conclude that ki% is unbounded at the 
poles. 

In a future work we will investigate the axially symmetric case, as initiated before [h5| 
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